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Abstract 

We study the effect of Debye screening of hypercharge when a net fermion number is 
reflected from a domain wall during a first order phase transition, which may be relevant 
for electroweak baryogenesis. We give a simple method for computing the effect of screen- 
ing within the diffusion approximation, whose results are compatible with those of a more 
elaborate treatment based on the Boltzmann equation. Our formalism takes into account 
the differences in mobility of different particle species. We believe it is conceptually simpler 
than other accounts of screening that have appeared in this context. Somewhat surprisingly, 
we find that Debye screening can actually enhance electroweak baryogenesis by a modest 
factor (~ 2). 



An important ingredient of the charge transport mechanism of electroweak baryogenesis 
|2j is the diffusion of net left-handed fermion number into the false vacuum, after quarks 
and leptons reflect from expanding domain walls during the electroweak phase transition. 
This left-handed asymmetry is what drives sphalerons to produce a baryon asymmetry It 
is desirable to understand how it propagates into the plasma because the more efficiently it 
does so, the longer it can drive baryon production before getting stopped by collisions and 
overtaken by the advancing wall. A number of papers have examined the diffusion problem 
from various angles @-f|. 

In the early days of the charge transport proposal, there were some worries that the 
chiral asymmetry, because it carries hypercharge, might be prevented from penetrating deep 
into the plasma due to Debye screening, since hypercharge is an unbroken gauged charge 
0. Indeed this would be the case if there were only one species of fermion carrying both 
chirality and hypercharge. It was argued however that when there are many species, the 



hypercharge can be screened without implying that chirality is also screened ||T0|, [TH| . In fact 
reference [[Uj concluded that Debye screening has no effect whatever on the rate of baryon 
production (up to small effects due to the difference between Fermi-Dirac and Bose-Einstein 
statistics. In the present work we reexamine the effects of screening and show that in fact 
it can have some effect on the baryon violation rate, actually making it somewhat larger 
than if one had neglected screening. The reason for the discrepancy in conclusions is that 
we allow for the various species to have different diffusion coefficients, which leads to some 
particles being more efficient screeners than others. We show how to explicitly compute 
the effects of screening by incoporating the gauge forces responsible for it into the diffusion 
approximation. In this formalism one can derive rather elegantly the final values of particle 
asymmetries after screening and diffusion, starting from some initial asymmetries produced 
at the domain wall. 

Our starting point, similarly to references || and |[LG|| , is the Boltzmann equation, in- 
cluding force terms for long range interactions. For each species of particle it is 

where i/i is the hypercharge of the ith species, E is the field strength and Cj is a collision term. 
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We take fa and rii below to refer to the net density of a given species, that is particle density 
minus antiparticle density. Integrating over momenta gives the usual continuity equation, 

Although the screening force appears to have been lost, it reappears in the diffusion equation 
when we make a derivative expansion of the current in terms of the density, because one 
must also add an Ohm's law term to J: 

Ji = -D t V ni + aiE, (3) 

where <7j is the conductivity and Di the diffusion constant of the ith species. 

There is a simple relation between the conductivity and the diffusion coefficient. The 
latter can be expressed in terms of the mean interaction time and the mean squared velocity 
(between collisions) as D = v 2 r/3, which for a relativistic particle is simply r/3. The 

— * 

conductivity on the other hand is the induced current divided by \E\, and the induced 
current is nv, where v is the average velocity due to the biasing caused by the electric field. 

— * 

After a collision the electric field causes the particle to change its momentum by p — yi\E\r 
in the time before the next collision, and v — p/p. The result is that 

n / n t2 / 0-33 bosons m 

a t = SnyMp = y lDl T j ^ (4) 

Now we want to imagine the situation where a macroscopically large domain wall, moving 

— * 

with a velocity v in the z direction, is the source of the particle densities. Thus Vn becomes 
dn/dz = n' . Since we are interested in solutions of the form n(z — vt), which would be static 
in the rest frame of the wall, we can replace dn/dt with — vn'. Furthermore the divergence 
of the electric field is given by Gauss's law, V • E — Y,Vi n i- So 

vn'i + Din" -aiJ2 Vj n j = °, ( 5 ) 

3 

or in vector notation on the space of species, 

vn' + Bn" -(y-n)a = 0, (6) 
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where D is the diagonal matrix with entries D{. This equation can be readily solved using 
the Laplace transform, 

N{s) = / dze- zs n{z). (7) 
Jo 

Making the definitions 

n = n(0); ft' = n'(0); 

f = vs + Bs 2 ; S = r 1 a; 

g = s"% + f _1 Dn , (8) 

equation ([]) becomes N — {y ■ N)Y, = g. But by dotting y into this equation one can solve 
for y ■ N and therefore for N, 

The solution for n(z) can then be obtained from the inverse Laplace transform 

n{z)= I^e zs N(s), (10) 

where the contour is taken to encircle the left half-plane (including an infinitesimal vertical 
strip of the right half-plane). To evaluate n(z) we need to identify the residues of the poles 
of (0) contained within the contour. Furthermore the initial values n and n' Q are related to 
each other by the physical requirement that the trivial constant solution does not exist and 
that the solution has the correct boundary condition vn(0) = n . So one must eliminate n' 
in favor of n to express the the solution in closed form. 

Let us start with the simplest situation, in which there is only one kind of particle. Then 

^-( 1 + 7^)(? + ¥)' 

which has poles at s = 0, s = -v/D and s± = -v/2D ± {v 2 + ADay) 1/2 /2D. But the 
residues at s = and s = Sq vanish, and s_ is the only one remaining within the contour. 
The resulting solution in z-space is n(z) = noexp(s-z). The effect of Debye screening is 
evident in this solution because the penetration depth is smaller than it would have been in 
the absence of the long-range force. In fact we can identify the Debye screening length as 

= (ay/D)- 1 ' 2 . (12) 
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When we take the limit of no screening, ay = 0, the solution of course reduces to that which 
has been previously found, n Q exp(— v z/D), showing that the penetration length gets longer 
as the wall gets slower, as one would expect from computing the time the wall takes to catch 
up with a diffusing particle. 

Next consider the case where there are many particle species, but they all have the same 
diffusion coefficient, D. Again N(s) has poles at s = 0, s = sq = —v/D and s = s±, 
but now s± is given by s± = —v/2D ± [v 2 + ADa ■ y) 1 ^ 2 /2D. In this case none of the 
residues automatically vanish; however the one at s = must be made to do so for physical 
consistency, since it would give a spatially constant solution. This gives the condition 

p.(n + °n j = 0, (13) 

where P (Q) is a projection operator onto the subspace orthogonal (parallel) to hypercharge, 

P.X-Q^X-lf). (14, 

Furthermore we get a relation between n and n' from the requirement that n(z) indeed 
gives n at z = 0, 

Q • n' = s_Q • n . (15) 
Thus the complete solution can be expressed in terms of no, 

n{z) = {e S0Z P + e s ~ z Q) ■ n , (16) 

which has a quite simple geometric interpretation: there is only screening in the direction of 
particle species space proportional to hypercharge. 

In the foregoing examples one could arrive at the solutions without recourse to the ma- 
chinery of Laplace transforming, but in the general case where there are several diffusion 
coefficients, the solutions would be harder to guess. We assume that there are M groups 
of particles classified by the values of their diffusion coefficients D a (for example different 
flavors of quarks have approximately equal diffusion constants because of their interactions 
with gluons). It is useful to define operators analogous to (III) that are nonzero only in the 
subspace of species with a given value D a : 

P^l„-Q«=(l„-4^), (17) 



where 



l a = lim ieCD-Da + te)- 1 ; 

e— >0 

a a = l a - a; y a = l a - y. (18) 

There are 2M physical solutions for the poles of N(s) which lie in the left-half plane. M of 
these are obvious generalizations of the previous unscreened solution, 

s a = ~v/D a . (19) 

The other M physical poles, which we shall call S{, occur at the negative roots of an (M+l)th 
order equation, 

sU(sD a + v)-^ff tt -if a '[l(sD b + v)=0 (20) 
The complete solution can be written as 

n{z) = }_^e a P a -n + 2^ °o. = — = Vb • n (21) 

a=l o,i,6=l ab ' y b 

where M ai is a square matrix with elements 

Mai = _ 2n l — (22) 
sfD a + Sit- 
in practice it is not no but rather the fluxes J(0) at the boundary which are known at the 
outset. However it is straightforward to prove the intuitive result that 

J(0) = vn (23) 



from eqs. (Q) and ( pT|) using the identities D a M ai Si = s { 1 — t>M ai , eq. fl2"0"D, and the electric 
field |£(0)| = J^y-ndz. 

It is trivial to verify that (pl|) satisfies the original differential equation; the nontrivial 
part is the determination of the roots Sj. If we assume that t> 2 /.D <C <? • y, which is the 
condition that the Debye screening length is much shorter than the usual diffusion length 
(and is always satisfied by the domain walls of the electroweak phase transition), it can be 
shown that two of the solutions to (|20|) are given by 

±(a a -y a /D a ) 1 / 2 + 0(v/D), (24) 
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of which the negative one is physical and gives rise to a single direction in the space of species 
which is screened. In the same regime of small velocities, the remaining M — 1 solutions are 
approximately the roots of the truncated version of eq. in which the first term (being 
of higher order in v) is omitted. This gives M — 1 values of Sj which are all of order v/D 
and therefore correspond to unscreened contributions. 

To illustrate the full solution we will consider a few special cases of interest. First let us 
assume that there are just M = 2 different diffusion coefficients, which could be useful for 
comparing the transport of quarks versus leptons, for example. Using eq. (|j) and ignoring 
the difference between the bosons and fermions, we find the simple result 

n{z)= (j2e~ Vz/Da Pa + ^{R2Vi - Rim) ® (yi/Ri - fa/R*)) -n Q , (25) 

\a=i v ' y ) 

where 

va-y v (yi-yi + y 2 ■ y 2 \ ^ 



Di°2 • m + £>2<?i -Vi y-y\D 2 D 1 
and 

Ra = (27) 

y ■ y 

We have omitted the screened contribution here, which is concentrated much closer to 
the bubble wall than the unscreened parts. It can be checked that this result reduces to our 



previous one ( |T6| ) in the limit that Di — ► D 2 . It is not correct in the limit that one of the 
D^s vanish, which corresponds to the decoupling of one of the particles, in the sense that 
it becomes immobile. This is because our assumption that v 2 /D <C o ■ y is then violated, 
but since we expect this condition to be satisfied in nature even for quarks, which have the 
smallest diffusion length, eq. (^B) is still the relevant one for considering the limit when one 
of the -D a 's is much smaller than the other one. 

The application motivating us to find this result was the case of differential reflection 
between chiral fermions and antifermions from the advancing domain wall as a means for 
producing the baryon asymmetry during the electroweak phase transition. In this mech- 
anism, the net production of baryon number is proportional to the integrated density of 
n b l -\- n l l , a linear combination of the left-handed baryon and lepton number densities. We 
can write this as 

POO 

AB = c dz u-n(z), (28) 
Jo 
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where u is a vector with components Ui = 1 for every left-handed particle and zero for all 
others, and c is some constant. It is easy to verify that u ■ y = 0, which stems from the fact 
that the baryon violation interaction giving rise to (^8|) is invariant under hypercharge gauge 
transformations. In the case of all equal diffusion coefficients Di, this has the consequence 
that AB oc u- hq because u is orthogonal to hypercharge, 

u ■ P ■ n = u ■ n . (29) 

This is the origin of the statement that hypercharge screening has no effect on the baryon 
asymmetry: we would have gotten the same answer if we had ignored screening altogether. 

However when the differences between diffusivity of particle species are taken into ac- 
count, it is no longer true that hypercharge screening has no effect because, for example, 

u ■ P a ■ n ^ u ■ l a • n , (30) 

because u is orthogonal to the full vector of hypercharges y, but generally not to its projec- 
tions y a onto the space of species with diffusion coefficient D a . 

Let us next examine the case where there are M = 3 distinct diffusion coefficients with 
the hierarchy Di D 2 < D 3 , which is the case for quarks and left-handed and right-handed 
leptons, respectively. Since Di is small we can expand in it. To zeroth order in Di, the 
integrated densities in front of the wall are given by 



f 

Jo 



. ' dzn(z) = (iD aPa+ ^:^f®M R >;f R >A .n , (31) 
'o Va=2 2/2 • 2/2/ D 3 + y 3 - y 3 / D 2 ) 

which is just the same result as (^): we recover the results for M = 2 when one of the 
particles becomes immobile (D\ — > 0). The 0(D{) corrections which involve the y\ direction 
either in the initial or the final density are 

D ADR 2 R 3 



£i|Pi+ ^ I D 2 (R 1 -1)/R 1 




( yi/Rl 
H2/R2 \ I - no, 



with 



(32) 

D = R 2 D 2 + R 3 D 3 - AD = D 3 — D 2 , (33) 
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and the _R a 's denned as in eq. From fl32|) we see that any particle densities induced 

by screening from an initial density of quarks (proportional to yi ■ n if we identify quarks 
with the 1-direction), as well as any quark density induced from initial densities of leptons 
(i/2,3 ' ^0) by screening, propagates only with the diffusion coefficient D\ of the quarks, which 
is much smaller than that of the leptons. Although this result may be intuitively obvious, 
the exact proportions by which these contributions are suppressed were not. 

To obtain more concrete results, we need to specify the values of the hypercharges and the 
diffusion coefficients. The space of species consists of 18 left-handed quarks, 9 right-handed 
up-type quarks, 9 right-handed down-type quarks, 6 left-handed leptons, 3 right-handed 
leptons, and let us say 2 Higgs doublets, counting colors and flavors. We will take the Higgs 
bosons to have the same diffusion coefficient as the left-handed leptons. The hypercharge 
vector can thus be written as 

V = {Vq, Vi\ Vr) 

( 1 4 2\ 
Vq = (18 x -; 9 x -; 9 x --J ; y q ■ y q = 22; 

yi = (6 x -1; 2 x -1) ; y r ft = 8; 

y r = (3x-2); y r ■ y r = 12. (34) 

Consider again the example with M = 2, where we assume that right and left handed leptons 
have equal diffusion coefficients (so that j7i ep ton = yi®Vr)- Then assuming that the inital flux 
consisted only of left handed bottom quarks with total density n q and right handed bottom 
quarks with density — n q and similarly for leptons with density ±n;, we would find 

AB = cv-\n q D q + mD{) - 3a;- 1 - - °« {D q n q - Dm) . (35) 

IL-Ug + lUJJl 

This result implies an enhancement in the flux due to the hypercharge screening by about 
30 per cent. It also illustrates the fact that regardless of screening, the particles with the 
largest diffusion coefficient are still the most efficacious for producing baryon number. 

The most interesting physical application is the situation where the initial state at the 
wall consisted only of left-handed tau leptons with density n T and an equal and opposite 
density — n T of right-handed tau leptons (since the process of reflection of fermions from the 
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wall conserves lepton number) |2|, |T^]. Neglecting the quark contribution entirely, we then 
find that 
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AB = cv " 1 T( 1 + rr(2A)/(3A.);-4 



~ -cv x Di 



(36) 



using that the ratio of diffusion constants of left- and right-handed leptons is approximately 
Di/D r = 1/4. The naive answer ignoring screening is AB = cDi/v. We thus see again the 
curious effect that screening increases the net amount of propagation of left-handed lepton 
number into the plasma. 
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Figure la Figure lb 

Figure 1. Figure la shows the hypercharge density in front of the wall. Figure lb shows the 
spatial distribution of the asymmetries of the left handed fermions. The solid line corresponds 
to the numerical solution of equation ( pi]) for the r^-asymmetry in the presence of screening 
and the short dash - long dash one to the r^-asymmetry in the absence of screening. The 
dashed line represents the eL,^L and //^-asymmetries {i = e, fi, r), which would be zero in 
the absence of the hypercharge screening. Distance is in units of T -1 . 



We show some results from the full numerical solution of the equation (|2T| ) corresponding 
to the above case of initial tau fluxes only, and with D q = 6/T, D\ = 110/T and D r = &D[. 
Figure la shows the hypercharge density in front of the wall. The width of the distribution 
is seen to correctly correspond to the inverse of the eigenvalue (^). Figure lb shows the 
distributions of the asymmetries of left handed leptons. Note in particular the distributions 
of the cl-i^l and neutrino-asymmetries that the hypercharge screening pulls out from the 
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plasma. It is the appearance of these new leptonic asymmetries, along with the fact that 
the Tl asymmetry is only slightly reduced by the screening, that explains the enhancement 
in the baryon production found above in the analytic formula ([36]). The precise numerical 
result for the final baryonic asymmetry is close to the analytic estimate: AB ~ l.Q9cv~ 1 Di. 

Let us finally caution the reader that there may be problems with the fundamental as- 
sumptions in the diffusion approach; Fick's law implicitly assumes that the reflected asym- 
metry current has a near-thermal distribution, whereas the initial flux emanating from the 
wall is rather strongly concentrated at low momenta |], 0, 0. It remains to be seen how 
much this will affect the above results of the diffusion computation. 

In conclusion, we have studied the diffusion of a reflected fermion current from the do- 
main wall during a first order phase transition, which is relevant for the charge transport 
mechanism of electroweak baryogenesis. In our treatment we included the effect of Debye 
screening of hypercharge and accounted for the differences in mobility of different particle 
species. Our main result is that in the case where the baryogenesis is driven by the flux 
of reflected tau leptons, Debye screening actually enhances electroweak baryogenesis by a 
modest factor ~ 1.7. 
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